QUANTIZATION OF THREE- WAVE EQUATIONS 
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Abstract. The subject of this paper is the consecutive procedure of discretization and quantization 
of two similar classical integrable systems in three-dimensional space-time: the standard three-wave 
equations and less known modified three-wave equations. The quantized systems in discrete space- 
time may be understood as the regularized integrable quantum field theories. Integrability of the 
theories, and in particular the quantum tetrahedron equations for vertex operators, follow from the 
quantum auxiliary linear problems. Principal object of the lattice field theories is the Heisenberg 
discrete time evolution operator constructed with the help of vertex operators. 



The tree-wave equations is the example of a completely integrable classical system in 2 -f 1 dimen- 
sional space-time. These equations apply to many physical systems. The aim of this paper is the 
formulation of an integrable quantum field theory corresponding to the three-wave equations. 

Given a classical theory with action S'[0], there is the universal prescription for the quantization. 
The amplitude between a state 4>in(j) at time ti and a state (j)out{ic) at time ^2 is 



'#>(i-.t2) = *o,.t(r) 

The perturbation-theory approach to the definition of D(/) we reject from every beginning. The practical 
way to define the measure in the Feynman integral is the discretization. It is to be understood as 



where the limit symbol stands for infinitely dense discretization of the space-time (r,t). 

Thus, the first step toward the field theory is the discretization of classical system. Only then 
we may look for a self-consistent Heisenberg quantum mechanics on the lattice and regard it as the 
regularized field theory. 

A schematic outlook of milestones of our method is the following. A classical theory is defined by 

the time dynamics for a field A, — A — f\A] (space-like degrees of freedom are omitted for brevity). 

at 

The dynamics is generated by a Hamiltonian, f[A] = {H, A}, where {, } are properly defined Poisson 
brackets. Corresponding discrete time dynamics is the evolution transformation A{t + At) — F[A{t)]. 
For brevity, we choose the scale At = 1. The key point is that making the discretization, we have 
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to take care of the Hamiltonian structure of the dynamics. The discrete time evolution must be a 
canonical transformation, it should preserve properly defined Poisson algebra of observable fields on 
space-like lattice. The quantum algebra of observables is a result of Dirac quantization of the Poisson 
algebra. The quantum evolution map A{t) — > A{t + 1) must be the an automorphism of the algebra 
of observables allowing one to define the Heisenberg evolution operator, 

(3) A{t + 1) = U A{t) . 

Such scheme is the algebraic realization of the discrete measure . Operator U acts in a representation 
space of the algebra of observables, the amplitude ([T|) in the Heisenberg form is 

» T-l T 

(4) {(|>out\U^\c|)^n)^ TT^'^* TT('^*l^l'^t-l) ■ 

The discretization of the space-time, construction of the quantum algebra of observables and the 
basis-invariant definition of Heisenberg evolution operator for the standard and modified three-wave 
systems is the subject of this paper. 

1. Three- WAVE equations 

We commence with a short reminding of the classical three-wave equations in three-dimensional 
space. In what follows, we use the short notations for the indices, 

(5) (j, J, k) — any permutation of (1, 2, 3) . 

1.1. Standard three- wave equations. Linear problem for the standard three- wave equations is the 
set of six differential relations 

(6) diTpj = A^jijji , z ^ j 

for three auxiliary fields V'i- Consistency of (O gives the equations for the six fields Aij-. 

(7) diAjk — AjiAik . 

These are the equations of motion for tree- wave resonant system (see e.g. [U H]), or the three- wave 
equations for the shortness. 

1.2. Modified three-wave equations. There is another type auxiliary linear problem [3^, the second 
order differential equations for the scalar auxiliary field cj): 

(8) {d^dj - A.jdj - AjA + A,,Aj,)(f> = , i^j . 
Consistency of ([5]) gives similar equations for the fields Aij , 

(9) d,Ajk = {A,j - Aik)iAjk - Aji) . 
We call them the modified three- wave equations. 
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1.3. Hamiltonians. It is convenient to use the single alphabetical indices instead of the numerical 
pairs. In this paper we will use the following convention for both systems: 

/-.QN ^12 — J ^13 — ' ^23 — ' 

A — A* A — A* 4oo — A* 

Note, our notations are not cyclic with respect to (1, 2, 3). 

Equations ([7]) and ([9| are extremum conditions for the action 

(11) S = j cPx {Ald:,A, - A*AA, + Ald^A, ~ V) 
where 

(12) V = A:A,AI - A.AIA, 
for the standard three- wave equations (l7|), and 

(13) V = {A^-A,){Al~A,){Al~Al) 

for the modified three-wave equations ([9]). For a moment, we ignore the reality conditions, the star in 
all these notations does not mean the complex conjugation. 

The time derivative dt and space derivatives dx^ dy for action pT|) may be chosen by 

(14) di=dt-dx, d2^-dt, d3 = dt-dy. 

Let r = {x,y) stands for the space-like vector. The standard Lagrange transform relating Lagrangians 
and Hamiltonians gives 

(15) n = Jd'r {AldyA, + AldxA, + V) , 
so that the equations of motion ([7]) and ^ are 

(16) 1^-^^'^^ = {^'^a(r,i)} , etc. 
where the Poisson brackets are defined by 

(17) \Al{v',t),A^{v,t)\ ^ 6{v' ~v) , v^a,b,c. 



Any other same-time bracket is zero. 

Relations pH]) and are just one of many possible conventions for the field notations and space- 
time separation. All such conventions in the continuous case are equivalent. The choice (jlOp and the 
signs in (jl4p at this moment may be considered as an odd decision of the author. 
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2. The Discretization 

Now we proceed to the discrete analogue of the standard and modified three- wave equations ()7|9|) . 
The straightforward discretization of three-dimensional space gives the cubic lattice 

(18) X = .TiGi + 2:262 + 0:363 , XiEM. n = 7iiei + 71262 + 71363 , Ui E Z 
The discrete analogue of derivative is the difference derivative, 

(19) a, ^ A, , A,(/.(n) =^ 0(n + e,) ~ ^{n) . 

The point is that we apply the straightforward discretization to the linear problems (j6l8p . discrete 
equations of motion must appear as the consistency conditions. 

2.1. Standard Thr66-wav6 6quations. The discrete linear problem corresponding to 

(20) AiVj,n = Ay(n)V'»,n , 

provides the discrete equations of motion 

^''^ + l-^.(n)A,-(n) ■ 

2.2. Modifi6d Thr66-wav6 6quations. The discrete linear problem corresponding to ([5]), 

(22) 0n+ei+ej " Qjj (n)(/'„+e, " Qii(n)0„+e, + (n)(3jj (n)(/)„ = , 

where 

(23) Qy(n) = 1 + Ay(n) , etc., 
provides 

, . „ ^„ Qji(n)g,fc(n) + Q,j-(n)Qjfc(n) - Q,j-(n)Qjj(n) 

(24) yjfe(n + ej) = ^ ■ 

3. Constant time surface and evolution 

Discrete equations of motion (|21l24p evidently define a sort of one-step evolution: the fields at point 
n + are expressed in the terms of the fields at point n. However, an explicit form of a space-like 
discrete surface and detailed definition of a discrete time corresponding to our choice (I14p needs some 
discussion. 

The key point is the geometrical structure of our discretization and the geometrical structure of the 
linear problems (I20|22|l . 

Let vector n in (fT8| stands for the vertex of the cubic lattice. The cubic lattice consists of vertices 
n, edges (n, n + e.;), and faces (n, n + e.;, n + 6j, n + + e^). The linear variable ipi.n of (pOI should 
be associated with (n, n + 6i)-edge, the linear variable 0n should be associated with n- vertex, and the 
fields Aij{n), Aji{n) should be associated with (n, n -I- 6^, n + 6j, n + e.; + ej)-face. This is justified by 
the structure of the linear equations, see Fig. [Hand Fig. [2l In what follows, symbol Av,n stands for 
the collection of fields on (v, n)-th face. 
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^l,n+e2 = "01,11 + ^2l(n)l/'2,n 
■02,n+ei = "02,11 + ^12(n)l/'l,n 



n + 62 4>l,n+e2 



Figure 1. Graphical representation of the elements of the linear equations (PO)) : lin- 
ear variables -01.11 are associated with the edges of cubic lattice, the fields Aa.n = 
(yli2(n), yl2i(n)) are associated with the face of cubic lattice. 



^ 0n+ei+e2 ~ Q210n+ei " Ql20n+e2 + Ql2Q210n — 



0n+e2 0n + ei+e2 



Figure 2 . Graphical representation of the elements of the linear equation : lin- 
ear variables (pn are associated with the vertex of cubic lattice, the fields Aa.n — 
((3i2(n), (52i(n)) are associated with the face of cubic lattice. 

Take now up the auxiliary linear problem (|20p and more detailed derivation of (|2ip . There are two 
ways to express V'l^n+ej+ea, '02,n+e3, 03,11+61+62 in the terms of -01,11, "02,11+61, '^a.n- The one way is 
to use the six relations 

A30i.„+e, = A3l(n + e2)03.n+62, A302.n = ^32(n)V'3,n, A201,„ = A2l(n)V'2,n, 

(25) 

Al-03,„+62 = ^13(n + e2)V'l,n+62, A2-03,n = A23(n)-02,n, Ai02^„ = A12 (n)-0i_„ . 

The second way is to use the six other relations 
(26) 

A3-02,n+6i = A32(n + ei)V'3,n+6i, A2VJi,„+63 = ^21 (n + e3)V'2,n+63 , A30i^„ = ^31 (n)03,n, 

A2V'3,n+6i A23(n + ei)?/'2,n+6i, Ai-02,n+63 = A12 (n + 63 )0i,n+e3 , Ai V'3,n = ^13 (n)01,n- 

Graphical representation for (|25p and (|26p is the left and right hand sides of Fig. [31 The collection of 
relations (pS)) and correspond to the faces of the cube (n, n + e^, n + + ej, n + ei + 62 + e-s) - 
the discrete space-time consistency is the consistency around the cube. 

The way to introduce the discrete time is to identify the discrete space-time fields in the left hand 
side of Fig. [3]with the discrete time t, and identify the discrete space-time fields in the right hand side 
of Fig. [3] with the discrete time t + 1. This exactly corresponds to p^ : 

(27) ei = et - e^r , 62 = -et , 63 = - , 
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"02,3 




Va,i2 



n + eg 




Vl,23 



Figure 3. Consistency around the cube for the linear problem of Fig[TJ Here -01 = 
"01,11, 01,2 = 0i,n+e2 7 "01,23 = V'i,n+e2+e3, etc. The left hand side corresponds to Eqs. 
((25)) . the right hand side corresponds to Eqs. ((26)) : Aa,n ~ (^I2(n),v42i(n)), etc., 
according to (fTO|) . 



and therefore 
(28) 



n = tet 



yey 



-ni, y ^ -ns , t = ni-n2 + n3 



The following table gives the correspondence between the initial Aij (n) - notations and the space-time 
notations, convention (1101) is taken into account: 



(29) 



^i2(n) = Aj2;,?;,i) = A, , A^^in) = Al{x,y,t) = A*^ , 

^I3(n + e2) = A^{x,y,t) = A^ , A^-^{n + e2) = Al{x,y,t) = Al , 
A3{n)^Mx,y,t)^A^, A^^ir^) = Al{x,y,t) = A* . 

Here the third columns are the shortened notations. The right hand side of Fig. [3] implies thus 

A,2{n + e-i) = A^ix,y - l,t + 1) =A^ , A^^in + e^) ^ Alix,y - l,t + 1) ^A* , 
(30) yli3(n) =ylb(x,2;,i + l) = ^ , A-^^in) = Al{x,y,t) ^A^* , 

A^sin + ei) = A^{x - 1, y, t + 1) = 3, , ^32(11 + ei) = A^x - l,y,t + 1) ^ A* , 

Here the third columns are the shortened notations as well. The consistency condition of (p5)) and 
([26]) . completely equivalent to Eq. ([2T|l . may be rewritten in shortened notations as 



A„ 



4 * 4- A"^ A 



A„ 



1 — A^|A^^ 



(31) <^ a:^ai{i-a,ai 



A* A* 



a:(i-a,^:)+a„a*(i-a,4*) 



A 



1 — ^f,^c 



, A^ — A^ + A*Ai, . 



In the absolutely similar way one may consider the discrete linear problem of Fig. [51 Corresponding 
picture is Fig. [H 
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023 



0123 



023 



0123 



Figure 4. Consistency around the cube for the hnear problem of Fig[21 Here cj) — 
01 = V'n+ei, 012 ~ 0n+ei+e2: 6tc. Convention (jlOp is taken into account, Aa.n - 
(Qi2(n),(52i(n)), etc. 

Introducing the same-time variables for the left hand side of Fig. [5] analogously to 

(5i2(n) = Ua(x,y,t) = Ua , Q2i(n) = Wa{x,y,t) = Wa , 

(32) Qi3{n + e2) ^Ub{x,y,t) ^Ub , Q3i{n + 62) ^ Wb{x,y,t) ^ Wb , 
Q23(n) = Uc{x, y, t) , Q32(n) = Wc(2:, t) ^Wc , 

and for right hand side of Fig. [?] analogously to ([5D1) . 

(3l2(n + 63) = Ua{x,y - 1, t + 1) = Ma , (52l(n + 63) = Wa{x, y - 1 , < + 1 ) = Wa , 

(33) (3i3(n) = Ufc(x,y,i + 1) = Ufc , (53i(n) = u;f,(a:, y, t + 1) = Wft , 

(323(n + ei) = Uc(a; - 1, y, < + 1) = Uc , Q32(n + ei) = Wc(a; - 1, y, < + 1) = Wc , 
we get the consistency condition in the shortened notations: 



(34) 



UaWb 


+ UbWc - 


UbWb 


Wa = 




WaWbUc 






Wc 


: 


UcWc - 


f WaWb - 


WaWc 


UaWa 


+ UbUc - 


UaUc 


Wb = 


UcWc - 


f WaWb - 


WaWc 






1 




Uc 






WaUbUc 




Wc = 


UaWb ^ 


f UbWc - 


UbWb 


UaWa 


+ UbUc - 


UaUc ' 









These relations are completely equivalent to Eq. (p4|) . 

The derivation of the equations of motion as the consistency of linear problem around the cube is 
similar to the consistency approach to the integrable equations on quad-graphs The consistency 
around the cube for two-dimensional quad-graph equations is an analogue of the Yang-Baxter equation 
- variables in the left and right hand sides of two hexagons of Fig. [3] are the same. The consistency in 
this paper implies the different variables in the left and right hand sides, it defines a map. Therefore 
it is an analogue of the local Yang-Baxter equation or the tetrahedral Zamolodchikov algebra [5l[6l[7], 
see next section. 
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Equations ([3T|l and (f34|) express the discrete space-time fields {Av{x, y,t+l) : v = a, b, c; x,y E Z} 
in the terms of the fields {Av{x,y,t) : v = a,b,c; x,y £ Z}. The space-like discrete surface is the 
collection of hexagons of Fig. [3] for all cc, y G Z and fixed t. This is the honeycomb lattice. Equations 
(pT|) and (jM]) therefore define the discrete time evolutions of the fields situated at the faces of the 
honeycomb lattices. 

Often in the literature the dual lattices are used. The fields are associated to the edges of dual 
three-dimensional lattice and to the vertices of its section - dual two-dimensional lattice. The lattice 
dual to the honeycomb one is called the kagome lattice. 

4. Linear problem as the zero curvature representation 

Equations ((5T|) and ([M)) are identically equivalent to Eqs. (^1]) and ((M)) correspondingly, but look 
more complicated since we reverse the time direction of 62 in Eq. (12711 . Our way to introduce time 
and space-like coordinates (|14l27p is not yet motivated. 

Equations (|3T|) define the map of variables ^„ = {A^,Al) to A^, — {A^,A*), v = a,b,c. Analogously, 
equation ([34]) define the map of variables Ay = w^) to Ay — {Uy,Wy), v — a,b,c. Denote both 
these maps by the symbol TZabc. Formally, TZabc is the operator acting in the space of functions of 

Aa ; Ab , Ac • 

(35) V$ = $(A,A,A) : (^7^abcO$^(A,A,A) $(X,A,X) . 

The point is that the maps ([?T|) and ([M)) are two basic set-theoretical solutions of the tetrahedron 
equation: 

(36) T^abcT^adeT^bdfT^cef = T^cefT^bdfT^adeT^abc ■ 

It may be verified straightforwardly, see the appendix. 

Take up now the question: why the maps (j3ip and p4p satisfy the functional tetrahedron equation. 
The reason is that the auxiliary linear problems defining the maps of dynamical variables from the 
left hand sides of Figs. [3] and |3] to the right hand sides are the correctly oriented zero curvature 
representations of three dimensional integrable model in discrete space-time. 

Let me demonstrate this statement for the map ([31]) . The linear equations for the single face, see 
Fig.m 

(37) V'l^n+es = "01, n + ^QV'2,n , V'2,n+ei = i^2,n + ^aV'l.n • 

may be rewritten in the matrix form as 

(38) f ) = X[A.] . f^^;"+-) , X[A.] = ( ] \ . 

The linear equations for all three faces of the left hand side of Fig. [3] may be written as 

(39) ('^^\= X[Aa] ( "^^'^1 , ( "^^'^ 1 = X[Ab] ( "^^'^^1 , f - X[Ac] ( "^^'^1 
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with the same matrix function X[A] (I38p . Iterating these matrix equations, we come to 



(40) 



( \ 

\ '03,12 / 



^12K]^13[A]^23[A 



^ -01,23 ^ 
^^2,3 



where are the three by three block-diagonal matrices, Xij coincides with ([38|) in (ij)-block 

and has the unity in complimentary block. For instance, 

^1 -Al \ 

(41) Xi2[Aa] - 

y 

Analogously, the right hand side of Fig. [3] provides 

( ^1 \ 



1 ^a^a 



(42) 



■02,1 
\^ "03,12 j 



X2z[Aa]XlMXi2[Ac 



' -01,23 ^ 
02,3 

v 



Thus, the consistency of the linear problem around the cube is the Korepanov equation [51 [B] 

(43) Xi2[Aa]XM]X23[Ac] = ^23pc]^13[A]^12[A] • 

Using definition (j35p . we may rewrite the Korepanov equation in the form similar to tetrahedral 
Zamolodchikov algebra [7j, 

(44) X,2[Aa]Xi3[Ab]X23[Ac] = nabcOX23[Ac]Xi3[Ab]Xl2[Aa]. 

Tetrahedron equation (|36|) is the equivalence of decompositions of the uniquely defined map Ay — Av 
(45) 

Xi2[Aa]X,3[Ab]X23[Ac]XiMd]X2Me]X3Mf]=X3Mf]X2Me]X23[Ac]X,4Ad]Xi4Ah^^^^ 

into two different sequences of elementary maps. 

The main difference between (|43)) and the local Yang-Baxter equation O [6] is that the numerical 
indices in (|43l44p correspond to the components of the tensor sum of one-dimensional vector spaces 
ipi. The Korepanov equation comes from the linear problem directly, therefore it is genuine multi- 
dimensional generalization of the Lax representation. 

There is no analogous matrix form for the zero curvature representation of the modified three- wave 
equations. One has to work directly with the sets of linear relations (j22|l . Nevertheless, the orientation 
of the faces in Fig. |4]is correct, and a set-of-linear-relations analysis similar to Eq. (j45| provides the 
set-theoretical proof of the corresponding tetrahedron equation [HI HU] . 

5. POISSON BRACKETS FOR THE FUNDAMENTAL MAPS 



The maps ([31]) and (|34|) define the discrete-time evolution on the honeycomb lattice. Evolution is 
the Hamiltonian one if it preserves Poisson brackets. 
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One may verify it for the map (|34p . if 



(46) 



and any other type bracket for Uv, Wy is zero, then 



Uv,Wv (■ = UvWv , V = a,b,c 



u-u, Wvf^ UyWv , V — a,b,c 



(47) 



and any other type bracket for Uy,Wv is zero. Therefore, ([34| is the canonical map |9]. Restoring 
the space structure according to p2l33p . we come to the whole system of same-time Poisson brackets 
conserved by the evolution: 



(48) 



Uy{x,y,t)wv{x,y,t) , 



any other type same-time bracket is zero. Poisson algebra with such structure of delta-symbols is 
called the ultra-local one. 

The Poisson brackets for the evolution ([3T|) are not ultra-local. To get the ultra-locality, we need 
to modify the discrete linear problem (|38p . Let there 



(49) 



X[A] 



K -A* 



, A^{A,A\K), K^^l-AA* 



\ A K ^ 

The modified map comes from the Korepanov equation (j43l 



A* - \k^AI + K^AIA^) , A^^K^ \KAa + K^A^Al) , 



(50) 



< - K, \KaA: + K,A,Al) , A, = K, \k,A, + K.AIA,) 



By definition P5)) . — 1 — A^A^ . Additional property of this map is 



(51) 



The map ([SO]) satisfies the functional tetrahedron equation and preserves the ultra-local Poisson brack- 
ets [TT] 



(52) 



-K,A,, 



any other type bracket is zero. 

The principal advantage of the Poisson structure is that it allows one to define the lattice actions with 
the help of functions generating the canonical transformations. This subject is technically complicated, 
we postpone it for future publications. 



11 



6. Quantization 

The Poisson algebra (j52|) is the quasi-classical limit of q-oscillator algebra [12] 

(53) AAt = 1 - g2N+2 ^ AtA = 1 - ^ = g^+iA , Atq^ = ^N-i^t . 

The Poisson bracket is the limit of commutator when — e^^ ^ 1, A ^ A, A^ — > A* and K . 

The Poisson algebra {u,w} — uw is the quasi-classical limit of the Weyl algebra [13] 

(54) uw = q'^wu 

with u u, w w when q — > 1. 

The whole algebra of observables, corresponding to the set of classical fields Av{x,y) is thus the 
tensor power of the local g-oscillator algebra A = (A, A''', g'^) for the quantized three- wave system; and 
the tensor power of the local Weyl algebra A — (u, w) for the modified three-wave system. After the 
quantization, the indices a, 6, c of the operators (or, more generally, the indices are (v, x, y), v = a,b, c, 
x,y G Z) stand for components of the tensor power. 

Quantum maps follow from the quantum linear problems. Quantized version of (j38p is [14j 



Here we replace the linear variables -04 by vectors [ipi) from a formal right module of the whole algebra 
of observables. Parameters Xa,fJ'a are C- valued spectral parameters, we introduce them for the sake 
of completeness. The consistency of linear problem of Fig. [3] (equivalent to the quantum Korepanov 
equation (|43)) ) may be solved with non-commutative coefficients, the answer is [11] [14] 



At„ = ^q-^^iq^^AK + ^g^"At,A,) , A, = ^q-^^iq^^A^ + -^g^»A,At,) , 

Afc q Ac AafJ-c 



(56) < 



At, = A,/Xcg^»+N^At, - At.At, , A, = ^gN.+N.A, _ A,A 



At, = ^,-N.(^N„At^ + h^q^.A,A\) , A, = l±q-^^iq^^Ac + -^q^^A^A,) , 
/i6 q fia AafJ-c 

Here q'^^'' = 1 - At;,Ab, in addition 

(57) ^N„+N. ^ ^N„+N, ^ ^N,+N. ^ ^N,+N. ^ 

One may verify, the map (|56p is the automorphism of the tensor cube of g-oscillator algebra (j53p . 
Quantized version of the linear equation ([22]) for the modified three-wave system is [10] 

(58) ^J4>n+ei+e2) ~ 9Wa|0n+ei) - VLal^ln+eo.) + UaVfa\4>n) = . 

Here Xa is a C-valued spectral parameter. Solution of the consistency condition, see Fig. [4] gives 



(59) 
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where 



>Ca >Cr _-Y ^a^c 



(60) A2 = — UftW^ + — u^w^ - g UftWj, , 

A3 ^ - grU^^W^ + ^f^U^ . 

One may verify, the map ([55]) is the automorphism of the tensor cube of the Weyl algebra 

Both automorphisms and (|59p satisfy the "functional" tetrahedron equation. In irreducible 
representations they are the internal automorphisms, 

(61) Tlal>c O $ = $ . 

Corresponding operators Robe satisfy the quantum (operator-valued) tetrahedron equations. Matrix 
elements of Rahc are functions of the spectral parameters Au, [i^ for (|56p and >f„ for ([59|). 

The local maps ([55)1 and ([55)1 define the evolution map on the honeycomb lattice via the identification 

(62) A = A,^,a(0 ^ A = A.a=,y(^ + 1) 

in accordance with (|29l30p and p2|33|) . The evolution map is the automorphism of the whole algebra 
of observables. Parameters A^,/iu for the g-oscillator model and «;„ for the Weyl algebra model should 
be (x, ?;)-independent. Then, in proper representations, the evolution is the internal automorphism 
given by an evolution operator, 

(63) $(i + 1) = C/$(i)[/"i . 

Matrix elements of U may be constructed with the help of matrix elements of local Robe- 

Examples of irreducible representations providing a "good" quantum mechanics are the following. 
For the g-oscillator algebra it is the case of real — e"'', < g < 1, and the Fock space: the Fock 
vacuum is defined by A|0} = N|0) = 0. The dagger of stands for the Hermitian conjugation, 
= N. If in addition — , — and in (|56p are unitary parameters, then Rahc and U are the 

well defined unitary operators. Matrix elements of Rabc are given in [11]. Evaluation (|4]) of in 
the framework of normal symbols gives the Feynman-type integral in the discrete space-time. In the 
quasi-classical limit h ^ Q [q^ = e^^), it may be shown 

(64) {A:^,\U^\AJ = j DATiA* e^^I^-^*! , 

where S'[^, A*] is the lattice action mentioned in the previous section, and the measure of integration 
is as well the lattice one. 

The proper quantum mechanical representation of the Weyl algebra is given by the modular dual- 
ization [15] . In addition to the given local Weyl pairs 

(65) u = , w = , \Q,P\ = \h => ^ e''^ , 
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it is necessary to consider the dual local pairs 

(66) u' = e^^ , w' = e^'S , q'^ = e"'^ . 

The sets of equations (j59p and similar equations for dual pairs define the kernel of Rabc unambiguously. 
It is known, in this case Rabc and U are well defined unitary operators. The PQ-symbol of in the 
quasi-classical limit fi, — > is 

(67) {Pout\U^\Q^n) ^ j 'DP'DQ e^^^P^'^^ 
where the measure and the action are the lattice ones. 



7. Conclusion 

The honeycomb lattice evolution map (|62|) defined with the help of local automorphisms ((56|) end 
(|59p . and the proper choice of the Hilbert space define unambiguously the evolution operator for 
the lattice approximation of corresponding quantum field theory. These quantum field theories are 
integrable since the underlying quantum auxiliary linear problems (|55l) and (j58[) provide the existence 
of the complete set of the integrals of motion [8l [161 HZl [HI [H] • 

The one-step evolution ([S^ is the discrete form of the Hamiltonian flow (115116^ . The quantum 
lattice Hamiltonian H is defined by 

(68) U = e- , 

where At = Ax — Ay is the lattice spacing parameter. In the continuous At ^ classical q ^ 1 limit 
the quantum lattice Hamiltonian becomes exactly (fTS]) . However, on the lattice At is finite (in this 
paper we used the scale At = 1), and therefore H is not a polynomial in the algebra of observables. 
The evolution operator is the principal object of the field theory rather than a Hamiltonian. 

The fundamental problem of the field theory is the calculation of the spectrum of evolution operator. 
This is the open question in three-dimensional models. Spectral equations for the evolution operators 
are not known yet (except for the two-dimensional limit of the three- wave system, [T9]). 

Let me conclude the paper with a discussion of the role of simplex equations in quantum field 
theory. The linear problem is the starting point of the integrability. The tetrahedron equation ([36]) 
and the Yang-Baxter equation for the two-dimensional models are the elementary consequences of 
linear problem as the zero curvature representation. Nevertheless, the simplex configurations may be 
considered as fragments of the space-time lattice, for instance 

(69) ^ ^ j ^'^'^ ^^2C^'?^3 ■ ■ ■ ^('/'l 02 I -Rl2 1 01 , 02 ) ("^l : '/'s I -^13 1 01 7 03 ) (<^2 . '/'s I -^23 1 02 > '/'s )/ ■ ■ 

i'0 eksm 

Here for brevity we consider a two-dimensional theory and the triangle configuration. Such partition 
function corresponds to the Feynman integral ([!]) with the discrete measure definition ^ and some 
particular discretization. The Yang-Baxter equation (and the ZJ-simplex equations in general) is the 
condition of Z-invariance: the partition function (j69p docs not depend on particular details of the 
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discretization, it is an invariant function of the boundary fields only. Thus the natural role of simplex 
equations in the lattice field theory is the conditions of self-consistent definition of the measure ^ . 
Acknowledgements I would like to thank all the participants of SIDE VII conference for fruitful 
discussions. 



References 

[1] V. E. Zakharov and S. V. Manakov, "Resonant interaction of wave packets in nonlinear media", JETP Lett. 18 
(1973) 243245. 

V. E. Zakharov and S. V. Manakov, "The theory of resonance interaction of wave packets in nonhnear media", Sov. 
Phys. JETP 42 (1975) 842-850 
[2] D. J. Kaup, "The inverse scattering solution for the full three dimensional three-wave resonant interaction", Physica 
ID (1980) 45-67 

[3] V. V. Mangazeev and S. M. Sergeev, "The continuous limit of the triple r-function model". Theoretical and Math- 
ematical Physics 129 (2001) 317-326 

[4] A. I. Bobenko and Yu. B. Suris, "Integrable noncommutative equations on quad-graphs. The consistency approach", 
Lett. Math. Phys. 61 (2002) 241-254 

V. E. Adler, A. I. Bobenko and Yu. B. Suris, "Classification of integrable equations on quad-graphs. The consistency 
approach." Comm. Math. Phys. 233 (2003) 513-543 
[5] J.-M. Maillet and F. Nijhoff, "Integrability for three dimensional lattice models", Phys. Lett. B224 389-396 
[6] R. M. Kashaev, I. G. Korepanov and S. M. Sergeev, "Functional Tetrahedron Equation", Theoretical and Mathe- 
matical Physics 117 (1998) 370-384 
[7] I. G. Korepanov, "Tetrahedral Zamolodchikov algebras corresponding to Baxter's L-operators" , J. Stat. Phys. 71 
(1993) 85-97 

[8] I. G. Korepanov, "Algebraic integrable dynamical systems, 2 4-1 dimensional models on wholly discrete space-time, 

and inhomogeneous models on 2-dimensional statistical physics" , Preprint |arXiv:solv-int/9506003| 
[9] S. M. Sergeev, "3D symplectic map", Phys. Lett. A 253 (1999) 145-150 
[10] S. M. Sergeev, "Quantum 2 + 1 evolution model", J. Phys. A: Math. Gen. 32 (1999) 5693-5714 
[11] V. V. Bazhanov and S. M. Sergeev, "Zamolodchikov's tetrahedron Equation and Hidden Structure of Quantum 

Groups", J. Phys. A: Math. Gen. 39 (2006) 3295-3310 
[12] E. V. Damaskinsky and P. P. Kulish, "Deformed oscillators and their applications", J. Soviet Math. 62 (1992) 
2963-2986 

[13] V. V. Bazhanov and N. Yu. Reshetikhin, "Remarks on the quantum dilogarithm" , J. Phys. A 28 (1995) 2217-2226 
[14] S. Sergeev, "Quantum curve in q-oscillator model", preprint nlin. 81/ 05 10048 , accepted in International Journal of 

Mathematics and Mathematical Sciences (2006) 
[15] L. D. Faddeev, "Discrete Heisenberg-Weyl group and modular group", Lett. Math. Phys. 34 (1995) 249-254 
[16] S. M. Sergeev, "Auxiliary transfer matrices for three-dimensional integrable models" , Theoretical and Mathematical 

Physics 124 (2000) 391-409 

[17] S. M. Sergeev, "Quantum matrices of the coefficients of a discrete linear problem", (Russian) Zap. Nauchn. Sem. 
S. -Petersburg. Otdel. Mat. Inst. Steklov. (POMI) 269 No. 16 (2000) "Vopr. Kvant. Teor. Polya i Stat. Fiz." 292- 
307, 370-371 Translation: "Coefficient Matrices of a Quantum Discrete Auxiliary Linear Problem", Journal of 
Mathematical Sciences 115(1) (2003) 2049-2057 
[18] S. Sergeev, "Integrability of g-oscillator lattice model", preprint Physics Letters A 357 (2006) 417-419 
[19] S. Sergeev, "Evolution operators for quantum chains". Preprint |arXivmlin.SI/0607055| 



15 



Appendix A. Verification of the functional tetrahedron equations 

The functional tetrahedron equation for the map (|3T|) may be verified with the help of Maple 9.5 
routine: 

restart ; 

R:=proc(a,b,c,var) # x stands for A, y stands for A~* 
local X,Y; 

X [2] : = X [b] * (1-x [a] *y [a] ) -x [a] *x [c] ; 

Y [2] : = y [b] * (1-x [c] *y [c] ) -y [a] *y [c] ; 

X [1] : = (x [a] * (1-x [c] *y [c] ) +x [b] *y [c] * (1-x [a] *y [a] ) ) / ( 1-X [2] *Y [2] ) ; 

Y [1] : =y [a] +y [b] *x [c] ; X [3] : =x [c] +y [a] *x [b] ; 

Y [3] : = (y [c] * ( 1-x [a] *y [a] ) +x [a] *y [b] * ( 1-x [c] *y [c] ) ) / ( 1-X [2] *Y [2] ) ; 

simplify (subs ([x [a] =X[1] ,y[a]=Y[l] ,x[b]=X[2] ,y[b]=Y[2] ,x[c]=X[3] ,y[c]=Y[3]] ,var)) 
end; 

TE:=var->R(l,2,3,R(l,4,5,R(2,4,6,R(3,5,6,var)))) 
-R(3,5,6,R(2,4,6,R(l,4,5,R(l,2,3,var)))) ; 

for k from 1 to 6 do TE(x[k]); TE(y[k]); od; 

Verification of the tetrahedron equation for the map is 
restart ; 

R : =pr oc ( a , b , c , var ) 
local U,W; 

U [1] : = (u [a] *w [b] +u [b] *w [c] -u [b] *w [b] ) /w [c] ; 

W [1] : =w [a] *w [b] *u [c] / (u [c] *w [c] +w [a] *w [b] -w [a] *w [c] ) ; 

U [2] : = (u [a] *w [a] +u [b] *u [c] -u [a] *u [c] ) /w [a] ; 

W [2] : = (u [c] *w [c] +w [a] *w [b] -w [a] *w [c] ) /u [c] ; 

U [3] : =w [a] *u [b] *u [c] / (u [a] *w [a] +u [b] *u [c] -u [a] *u [c] ) ; 

W [3] : = (u [a] *w [b] +u [b] *w [c] -u [b] *w [b] ) /u [a] ; 
simplify (subs ( [u [a] =U[1] ,w[a]=W[l] ,u[b]=U[2] ,w[b]=W[2] ,u[c]=U[3] ,w[c]=W[3]] ,var)) 
end; 

TE:=var->R(l,2,3,R(l,4,5,R(2,4,6,R(3,5,6,var)))) 
-R(3,5,6,R(2,4,6,R(l,4,5,R(l,2,3,var)))) ; 

for k from 1 to 6 do TE(u[k]); TE(w[k]); od; 
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